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\  Introduction 

In  the  last  few  years,  several  investigators  WTttvW-l- _ 

have  studied  the  idea  of  estimating  lung  volume  by  measuring  the  dimen¬ 
sional  changes  of  the  chest  and  abdomen  during  respiration.  The 
prospect  of  being  able  to  know  the  value  of  lung  volume  from  information 
obtained  noninvasively  was  intriguing. 

The  most  common  method  of  measuring  the  dimensional  changes  is  to 

use  magnetometer  pairs.. _ A  more  complete  description  of  the  principles 

of  operation  can  be  found  elsewhere  [7].  -^For  the  purpose  of  this 
report,  it  is  sufficient  to  note  that  magnetometer  pairs  generate  a  vol¬ 
tage  that  is  proportional  to  the  change  in  their  separation  distance. 

The  relationship  between  voltage  and  the  change  in  separation  distance 
is  essentially  linear  over  the  separation  distances  measured  in  this 
study. ( 

If  one  is  to  adequately  infer  lung  volumes  from  the  measurement  of 
dimensional  changes  of  the  thorax  and  or  abdomen,  the  following  ques¬ 
tions  must  be  addressed: 

1) .  What  dimensional  changes  should  be  measured? 

2) .  How  are  the  dimensional  changes  related  to  lung  volume 


changes  ? 

3) .  How  reliable  can  the  dimensional  changes  be  measured? 

4) .  Does  the  relationship  between  dimensional  changes  and 

lung  volume  change  from  breath  to  breath,  or  with  the 
respiratory  maneuver? 


5).  Does  Che  relationship  between  lung  volumes  and 

dimensional  changes  differ  from  individual  to  individual? 

We  will  attempt  to  answer  many  of  the  above  questions,  and  will 
speculate  on  the  answers  to  the  others.  Before  describing  our  current 
experiments,  however,  it  is  worthwhile  to  review  recent  measurements 
with  magnetometers. 

The  use  of  magnetometers  in  studying  pulmonary  ventilation  was  pro¬ 
posed  by  Mead  et  al.  [7],  In  this  pioneering  work,  Mead  and  co-workers 
built  upon  an  earlier  study  [6]  in  which  they  proposed  that  the  abdomino¬ 
thoracic  cavity  could  be  treated  as  a  two  degree  of  freedom  system  (rib 
cage  and  abdomen),  and  that  for  isovolume  maneuvers,  volume  change  is 
nearly  linearly  related  to  changes  in  anteroposterior  (AP)  diameters. 
They  showed  [7]  that,  after  calibration,  the  sum  of  magnetometer 
measured  AP  rib  cage  and  abdominal  diameters  reflect  lung  changes 
closely.  They  recommended,  therefore,  that  magnetometer  measurements  of 
minute  ventilation  would  be  useful  where  conventional  (eg.  spirometric) 
techniques  were  inconvenient. 

Gilbert  et  al.  [3]  used  the  method  proposed  by  Konno  and  Mead  to 
investigate  breathing  patterns  during  CO2  inhalation.  The  major  change 
was  their  introduction  of  a  new  "calibration  procedure"  whereby  the  two 
magnetometer  signals  (chest  and  abdomen)  were  "scaled"  graphically  by- 
superimposing  two  breaths  (an  abdominal  breath  and  a  chest  breath).  In 
this  way,  they  found  that  the  "scaled-summed"  magnetometer  signal  was 
linearly  related  to  the  spirometer  volumes.  The  magnetometer-measured 
data  was  then  used  to  construct  tidal  volume  -  ventilation  curves  for 
ten  subjects . 
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In  a  later  study,  Gilbert  et  al.  [4]  used  magnetometer-measured 
tidal  volumes  to  show  that  conventional  methods  utilizing  noseclips  and 
mouthpieces  alter  the  pulmonary  parameters  of  :  respiratory  frequency, 
ventilation  and  tidal  volume. 

Grassino  and  Anthonisen  [5]  used  magnetometers  to  examine  the  degree 
of  distortion  of  the  chest  wall  at  functional  residual  capacity  (FRC) 
during  both  high  resistive  inspiration  and  under  external  lateral 
compressions.  In  addition,  they  used  magnetometer  pairs  to  describe 
thoracic  shape  alterations  while  simultaneous  regional  volume  distribu¬ 
tions  were  measured  with  a  Xeon  technique.  It  is  important  to  note  that 
their  studies  were  for  isovolume  maneuvers  and  that  the  results  of  Ronno 
and  Mead  [6]  were  applicable. 

Ashutosh  et  al.  [2]  used  magnetometer-measured  AP  diametrical 
changes  to  study  breathing  patterns  in  both  normal  and  COPD  (Chronic 
Obstructive  Pulmonary  Disease)  patients.  They  showed  qualitatively  that 
both  abdominal  and  chest  motions  where  synchronous  with  spirome  t r i  ca  1 1 y 
measured  breathing  in  all  10  normal  subjects  and  in  7  of  17  COPD  pa¬ 
tients.  In  the  other  10  COPD  patients,  the  chest  motion  was  found  to  be 
synchronous  with  spirometric  volume,  but  the  abdominal  signal  was  asyn¬ 
chronous.  It  was  further  shown  that,  in  general,  the  asynchronous  pat¬ 
tern  corresponded  to  a  poorer  patient  prognosis.  They  concluded  that 
recognition  of  this  type  of  breathing  pattern  could  be  most  helpful  in 
initial  patient  assessment.  Also,  in  an  interesting  application,  they 
used  the  magnetometer  signals  to  instruct  patients  to  improve  their 
breathing  patterns  by  matching  the  magnetometer- measured  breathing  pat¬ 
terns  with  normal  (desired)  pattterns.  This  method  was  used  to  help 
wean  patients  off  of  ventilatory  assist  devices. 
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Sharp  et  al.  [10],  using  what  they  call  the  "Konno  and  Mead  method 
of  thoracoabdominal  partitioning  of  breathing",  looked  at  81  normal 
subjects  to  investigate  whether  variations  related  to  sex  and/or  age 
differences  exist.  Using  two  magnetometer  pairs  (chest  and  abdomen), 
they  found  no  major  differences  in  relative  contributions  of  rib 
cage  /abdominal  breathing  between  men  and  women,  or  between  young  and  old 
during  any  respiratory  act.  In  addition,  they  pointed  out  two  most  im¬ 
portant  points:  first,  that  for  rapid  ventilatory  maneuvers,  the  approxi¬ 
mately  linear  relationship  between  volume  change  and  AP  diameter  no 
longer  prevails,  although  preliminary  studies  suggested  that  in  these 
ranges  magnetometer  based  information  is  still  qualitatively  useful ;  and 
second,  that  phase  lags  in  the  lateral  and  AP  diameter  changes  (rib  cage 
and  abdomen)  render  useless  any  attempts  to  interpret  phase  relation¬ 
ships  during  maximal  voluntary  ventilation. 

In  summary,  we  notice  that  the  aforementioned  investigations  have 
used  two  magnetometer  pairs  to  qualitatively  investigate  the  roles  of 
the  rib  cage  and  abdomen  in  breathing.  Now,  we  mention  several  investi¬ 
gators  who  report  quantitative  results  on  inferring  lung  volumes  from 
magnetometer  measured  diametrical  body  surface  changes. 

Stagg  et  al.  [11],  again  using  two  magnetometers,  introduced  both  a 
new  calibration  method  and  a  volume  model  from  which  tidal  volumes  could 
be  inferred.  They  showed  that  it  is  possible  to  calibrate  the  magneto¬ 
meters  accurately  during  spontaneous  breathing.  However,  they  (like 
Sharp)  pointed  out  that  there  was  no  reported  evidence  stating  that  chest 
wall  displacements  are  linearly  related  to  volume  at  the  extremes  of 
vital  capacity.  Thus,  they  concluded,  that  magnetometer  measurements 


should  be  used  within  "moderate"  volume  ranges.  They  suggest  the  techni¬ 
que  as  an  accurate  means  of  measuring  tidal  volume  and  the  time  compo¬ 
nents  of  individual  breaths. 

They  suggest  four  possible  sources  of  error  in  the  magnetometer 
based  method:  first,  the  calibration  procedure;  second,  accurately  de¬ 
fining  the  change  in  respiratory  phase  (expiration  to  inspiration); 
third,  compression  and  decompression  of  thoracic  gas  at  higher  rates  of 
ventilation,  and  fourth,  cases  of  abnormal  abdominothoracic  distortion. 

Robertson  et  al.  [9]  postulated  three  volume  models  (the  first  of 
which  was  analagous  to  Stagg)  to  quantitatively  evaluate  lung  volume. 
They  propsed  that  two  additional  magnetometer  pairs  (placed  laterally  at 
the  same  level  as  the  AP)  be  used  in  the  model.  They  showed  that  a  four- 
magnetometer  elliptical  cylinder  model  gave  the  best  results  in  quiet 
breathing  and  in  vital  capacity.  They  also  point  out  that  at  the 
extremes  of  lung  volume,  the  method  may  break  down.  Robertson  reports 
an  everpresent  counterclockwise  "looping"  of  the  estimated  volume  at  all 
tidal  volume  ranges.  They  suggest  two  possible  reasons:  first,  "that 
different  levels  of  the  chest  may  behave  differently  in  relation  to  the 
magnetometer  between  inspiration  and  expiration;  they  may  lag  behind  or 
precede  movements  of  the  magnetometers”.  Second,  the  shift  of  blood  to 
and  from  the  extremities  and  thorax  may  be  a  factor.  They  conclude, 
however  ,  that  this  method  may  be  particularly  well  suited  to  studies  of 
respiratory  control  and  patient  monitoring. 

Ackerman  [1],  using  the  method  of  Robertson,  automated  the  volume 
measurement  on-line.  He  reports  that  breaths  per  minute,  average  tidal 
volume,  and  minute  volume  can  be  inferred  and  displayed  at  15  second 
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intervals  while  monitoring  a  patient.  He  reports  results  are  accurate  to 
within  10%  of  spirometric  techniques. 

Vawter  [12],  tested  seven  volume  models  (one  being  the  same  as  the 
Robertson  model)  and  showed  that  any  two  or  three  parameter  model  is 
adequate  to  predict  volume  for  the  breath  from  which  the  model  constants 
are  determined.  He  further  reported  that  calibration  should  be  done 
with  a  complex  respiratory  maneuver  if  the  model  is  required  to  predict 
such  complex  maneuvers. 

Although  Vawter  suggests  that  no  one  dimensional  measurement  is 
sufficient  to  predict  volume  changes,  he  does  report  that  two  anatomical 
sites  (AP  chest  and  AP  abdomen)  individually  correlated  85%  or  higher 
with  spirometric  volume.  He,  as  did  Robertson,  noted  the  "looping"  or 
hysteresis  of  respiratory  movement  and  suggested  the  need  to  consider  it 
in  any  further  studies . 

Melissinos  et  al.  [8]  studied  changes  in  abdominothoracic  shape 
during  forced  vital  capacity  (FVC)  maneuvers.  Using  four  magnetometers 
at  different  sites  than  aforementioned,  they  reported  that  at  the  AP 
xiphi-sternal  junction  and  AP  manubrium  site  that  diametrical  changes 
with  volume  are  useful  indices  of  the  motion  of  the  anterior  chest 
during  FVC.  They  also  note  that  measurements  are  accurate  (from  iso¬ 
volume  calibration  [6])  during  spontaneous  breathing  and  slow  respira¬ 
tory  manuevers  (20-80%  vital  capacity).  They  found  changes  in  lateral 
xiphi-sternal  magnetometers  quite  variable  with  subject.  Also  they 
report  that  AP  abdominal  changes  may  not  be  representative  of  the  over¬ 
all  movement  of  the  anterior  abdominal  wall.  Finally,  they  demonstrated 
substantial  nonuni formites  in  regional  abdominothoracic  dimension  changes 
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during  FVC,  and  that  overall  chest  wall  volume  displacement  cannot  be 
accurately  represented  by  the  two  common  magnetometer  positions  (rib 
cage  and  abdominal)  during  FVC. 

METHODS 

In  order  to  answer  the  questions  posed  above,  we  have  conducted  a 
series  of  experiments  in  which  we  measured  simultaneously  both  dimen¬ 
sional  changes  and  changes  in  lung  volume. 

Four  magnetometer  pairs  were  used  to  measure  diametrical  changes  of 
the  abdominothoracic  cavity  at  eight  anatomical  sites.  Consistent  with 
previous  studies,  the  midline,  both  AP  and  laterally,  was  chosen  as  a 
locus  of  possible  placement  sites.  We  note  that  the  lateral  magneto¬ 
meters  were  placed  just  anterior  to  the  the  latissimus  dorsi  muscle. 

The  eight  sites  chosen  were  :  AP;  ,  just  superior  to  the  sternal  arch; 

M2,  at  the  level  of  the  xiphoid;  M3,  midway  between  the  xiphoid  and 
umbilicus,  and  M4 ,  just  inferior  to  the  umbilicus.  Laterally:  M5,  at 
the  level  of  the  fourth  rib;  Mg,  just  inferior  to  the  xiphoid  level;  My, 
midway  between  Mg  and  Mg,  and  Mg  at  the  level  of  the  umbilicus  (Fig.l). 

Eleven  normal  subjects  (ages  19-29)  with  no  prior  pulmonary  func¬ 
tion  testing  experience  were  studied  (table  1).  The  magnetometer  pairs 
were  taped  securely  in  the  above  mentioned  post  ions.  Care  was  taken 

i 

I 

that  the  long  axes  (yj  and  y£  see  Fig.  2)  were  parallel  to  avoid  i 

| 

rotational  effects  [12].  Proper  alignment  was  obtained  at  the  postion 
which  generated  a  global  minimum  in  voltage  when  one  magnetometer  was 
rotated  with  respect  to  the  other. 

Standing  erect,  each  subject  performed  two  separate  breathing  maneu¬ 
vers:  "quiet  breathing"  and  "forced  breathing"  (ie.,  one  vital  capacity 
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maneuver  was  performed  at  the  middle  of  a  quiet  breathing  sequence).  A 
total  of  eight  breathing  tests,  each  of  forty-five  seconds  duration  were 
performed  by  each  subject.  Thus,  two  "quiet"  and  two"forced"  breathing 
tests  were  performed  at  each  of  the  eight  magnetometer  placement  sites. 

A  given  test  consisted  of  the  subject  performing  the  particular 
maneuver  by  breathing  into  a  spirometer  (Model  840,  Ohio  Medical 
Products)  with  his  nose  clamped.  The  subject  was  instructed  to  minimize 
all  unnecessary  body  motions.  The  spirometer  and  the  four  magnetometer 
signals  were  recorded  simultaneously  as  voltage  (output)  versus  time 
(Fig.  3).  A  Mine  11  Computer  (Digital  Equipment  Company)  was  used  to 
convert  the  five  analog  signals  to  digitized  form,  and  then  to  store  the 
data  on  floppy  disk  for  later  analysis.  All  programs  were  written  in 
BASIC  (Appendix  1),  and  the  data  sampling  rate  was  6  6.67  samples  per 
second  (ie.  13.33  samples  per  channel  per  second). 

Quantitative  data  analysis  consisted  of  two  approaches:  First,  in 
order  to  investigate  the  relationship  between  spirometric  data  and  the 
diametrical  changes  of  the  chest  and  the  abdomen,  as  well  as  the  cross 
relationship  between  various  anatomical  sites  (Mj-Mg),  correlation  coef¬ 
ficients  (0 )  for  the  respective  data  were  calculated.  Second,  fourier 
analysis  was  used  to  represent  the  data  (ie.  a  given  breath)  as  a  sum  of 
sinusoidal  components  to  study  the  influence  of  "phase"  and  "amplitude" 
on  the  correlation.  Two  authors  [8,9]  have  speculated  that  "phase  dif¬ 
ferences"  may  influence  how  well  magnetometer  based  techniques  can  be 
used  to  infer  volume  information.  In  addition,  the  spectral  analysis 
allowed  us  to  ascertain  whether  the  signals  could  be  modeled  by  a  simple 
harmonic  function. 
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Since  dimensional  changes  during  expiration  and  inspiration  may  be 
different  [9,12],  each  portion  of  the  breathing  maneuver  was  analyzed 
separately. 

RESULTS 

The  correlation  coefficients  relating  the  spirometric  results  to 
each  of  the  eight  magnetometer  placement  positions  are  given  for  each 
subject  (Appendix  2).  It  can  be  seen  from  these  tables  that  for  each 
subject,  at  least  one  AP  position  had  a  correlation  with  the  spirometric 
data  of  greater  than  0.9.  In  all  subjects,  except  subject  #9,  the  M} 
location  has  the  highest  correlation  coefficient.  Subject  #9  is  an 
abdominal  breather  which  illustrates  the  need  to  allow  for  different 
subject  types.  Tables  2,  3,  4  and  5  list  the  averaged  correlation 
coefficients  +  standard  deviations  for  correlation  between  spirometer- 
magnetometer  position  and  cross  correlations  between  the  magnetometer 
positions.  Gross  observations  from  tables  2,  3,  4  and  3  show  that 
the  AP  positions  in  each  case  correlate  higher  than  their  lateral 
counter  parts.  As  can  be  seen,  the  Mj ,  position  tends  to  correlate  very 
well  (>0.95)  with  the  spirometer.  It  is  important  to  note  that  this 
high  correlation  is  coupled  with  a  very  low  standard  deviation.  On  the 
other  hand,  we  note  poor  correlation  between  the  spirometer  and  the  M2 
and  Mg  positions  (ie.  AP  and  lateral  xiphoid  levels).  Not  only  are  the 
correlations  (.1610  ^<.4638)  very  poor,  but  the  scatter  of  the  data  is 
reflected  by  the  large  standard  deviations  (.3818  -  .5982).  Also,  it  is 
seen  that  negative  correlations  arise  in  the  xiphoid  data  (subjects  6,7 
in  the  AP  position  and  subjects  6,9,10,11  in  the  lateral  position). 
Determination  of  the  cause  of  this  poor  correlation  would  require 


9 


alternate  measurement  of  the  dimensional  changes.  It  could  be  that 
paradoxical  breathing  is  the  cause  or  that  the  magnetometers  in  this 
position  were  subjected  to  rotational  motions,  to  which  the 
magnetometers  are  sensitive  [12]. 

The  only  cross  correlation  between  magnetometer  sites  (AP  and 
lateral  considered  separately)  with  a  80%  or  higher  correlation  is 
between  the  M3  and  M4  sites  (84%  and  83%).  Also,  note  that  the  corre¬ 
lation  between  magnetometer  pairs  is  much  lower  than  between  the 
magnetometers  and  the  spirometer.  If  the  correlations  were  high, then 
the  measurements  would  be  redundant  and  one  of  the  pairs  could  be  elimi¬ 
nated. 

The  correlation  coefficients  reported  above  were  calculated  from 
data  sets  of  seven  to  twelve  breaths  (45  seconds)  for  each  of  the  two 
maneuvers.  In  addition,  we  looked  at  correlation  coefficients  for  indi¬ 
vidual  breaths,  and  at  the  inspiratory  and  expiratory  portions  of  single 
breaths.  We  found  that,  in  general,  the  correlation  coefficient  for  an 
individual  breath  was  higher  than  that  of  the  total  breathing  sequence, 
and  as  the  number  of  breaths  increased  toward  the  total  for  that  indivi¬ 
dual  test,  the  correlation  coefficients  approached  that  of  the  total. 
These  results  are  to  be  expected  if  there  are  breath  to  breath  varia¬ 
tions  in  the  signals  (and  if  the  signals  have  only  a  small  component  of 
random  noise).  No  significant  patterns  were  observed  in  this  result 
(Table  6). 

In  table  7  we  show  the  correlation  coefficients  for  each  subject 
(AP  positions  only)  for  six  individual  breaths  from  a  "quiet  breathing" 
test.  As  can  be  seen  from  the  averages  +  standard  deviations,  position 
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Mj  has  a  high  average  correlation  (50.93)  and  a  low  standard  deviation 
(.002  -  .055).  As  in  the  above,  the  abdominal  breather  (subject  #9)  had 
the  lowest  M}  correlation.  Also,  the  M2  position  had  the  greatest  varia¬ 
bility  and  typically  the  largest  scatter  in  the  date.  We  observe, 
therefore,  that  breath  to  breath  variations  do  not  appear  significant 
except  in  the  M2  position. 

In  the  "forced  breathing"  manuever,  the  subject  was  instructed,  at 
a  particular  instant,  to  (on  his/her  end  resting  expiratory  volume) 
inhale  as  deeply  as  possible  (i.e.  to  maximum  inspiratory  level)  and 
then  exhale  totally.  This  can  be  seen  in  figure  3  in  the  eighth  breath. 
In  each  of  the  above  cases,  the  correlations  were  calculated  based  on 
either  averaged  breathing  or  on  quiet  breaths.  In  table  8,  we  show  for 
each  subject  the  AP  correlation  coefficients  calculated  from  a  single 
"forced  breath".  As  can  be  seen,  observation  of  the  forced  breathing 
results  showed  no  remarkable  differences  from  those  of  quiet  breathing. 

To  investigate  the  dependence  of  volume  on  magnetometer  measured 
dimensional  changes  we  plotted  one  versus  the  other.  In  most  every 
case,  hysteresis  was  present  showing  differences  between  expiration  and 
inspiration.  We  generated  these  plots  (spirometer  vs  M£,  i  =  1,2. ..8) 
for  three  breaths  for  each  subject.  In  general,  of  the  AP  positions,  it 
was  the  M^  position  which  gave  the  least  hysteresis  and  the  M2  position 
which  corresponded  to  the  most  hysteresis  (see  Fig.  4).  We  note, 
however,  that  considerable  variation  in  hysteresis  was  seen  on  a  breath 
to  breath  basis  (see  Fig.  5).  Similar  results  were  observed  laterally 
(Fig.  6).  The  degree  of  hysteresis  is  certainly  reflected  in  the  calcu¬ 
lated  correlation  coefficients:  this  can  be  seen  in  figure  7  where  three 
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different  plots  are  given  with  the  corresponding  correlation  for  that 
breath.  Thus,  we  felt  it  beneficial  to  separate  the  breath  into  its 
inspiratory  and  expiratory  portions  and  calculate  the  respective 
correlation  coefficients.  We  show  the  results  calculated  over  a  single 
quiet  breath  in  table  9.  It  is  worthwhile  to  note  that  in  each  case 
where  the  total  breath  correlation  is  0.99  or  greater,  the  inspiratory 
and  expiratory  results  are  essentially  the  same,  as  expected.  In  the 
case  shown,  it  is  interesting  that  both  subjects  for  which  higher 
abdominal  correlations  occurred,  (subject  #9  and  subject  #11  for  this 
breath)  the  expiratory  results  were  correlated  higher  than  the 
inspiratory  (for  all  >  0.99).  Conversely,  the  other  nine  subjects 
tended  to  have  higher  inspiratory  results  in  cases  where  the  results  did 
differ . 

The  fourier  analysis  (Appendix  3)  reveals  that  the  amplitude  of  the 
first  harmonic  was  normally  7  to  20  times  as  large  as  those  for  higher 
harmonics  for  both  the  spirometer  and  the  Mi  magnetometer.  The  other 
magnetometers  would  not  be  adequately  described  using  a  single  harmonic. 
In  figure  8  we  show  a  M2  plot  for  a  single  breath,  and  the  need  for  seven 
harmonics  to  adequately  describe  the  curve.  In  fact,  often  for  the 
xiphoid  level,  the  first  harmonic  was  not  dominant. 

We  also  give  in  table  10  the  calculated  correlation  coefficients 
for  the  breath  for  which  the  fourier  results  are  given.  Generally,  it 
is  seen  that  in  cases  of  high  correlation  the  phase  difference  between 
spirometer  and  magnetometer  position  is  lower.  Looking  at  the  magneto¬ 
meter  position  which  had  the  most  dominant  first  harmonic  (which  tended 
to  coincide  with  the  highest  correlation  for  that  subject)  we  see  no 
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consistent  pattern  as  to  the  dimensional  changes  lagging  or  leading  the 
spirometric  results.  Thus,  we  can  only  say  that  variations  from  subject 
to  subject  are  observable  in  this  regard. 

DISCUCSSION 

Unfortunately,  to  date  there  exists  no  literature  to  which  we  can 
directly  compare  our  results.  There  are,  however,  reported  findings 
and/or  speculations  to  which  our  results  may  be  addressed. 

Robertson  et  al.  [9]  and  Melissinos  et  al.  [8]  suggest  that  dif¬ 
ferent  levels  of  the  chest  and  abdomen,  respectively,  behave  differently 
in  relation  to  the  magnetometer  -  volume  results.  We  point  out  that 
each  of  the  eight  anatomical  sites  considered  herein  did  yield  different 
results  in  correlation,  hysteresis  and  phase  relationships.  Thus  we 
have  shown  that  although  within  a  given  anatomical  region  correlations 
may  be  similar  (eg.  M3  and  M4  vs  S  in  tables  2,3),  each  site  considered 
did  behave  differently.  Coupling  our  results  with  those  reported  by 
Vawter  [12]  further  substantiates  this  point. 

Sharp  et  al.  [10]  and  Robertson  et  al.  [9]  both  mention  that  the 
difference  in  phase  between  volume  and  dimensional  changes  may  influence 
the  usefulness  of  the  magnetometer- found  data.  Indeed  this  may  well  be 
true  but  may  be  difficult  to  quantify.  Only  the  Mj  position  yielded 
results  in  which  a  dominant  first  harmonic  was  seen,  and  as  pointed  out 
above,  subject  to  subject  variability  renders  this  analysis  quite  use¬ 
less  in  generating  a  general  conclusion. 

We  conclude  by  emphasizing  comments  by  Robertson  et  al.  [9]  and 
Ashutosh  et  al.  [2]  in  that  the  usefulness  of  magnetometers  in  studying 
pulmonary  function  may  indeed  lie  in  the  realm  of  patient  assessment  and 
patient  monitoring.  Certainly,  information  from  and  M3  or  M4 


positions  is  qualitatively  useful.  However,  useful  quantitative 
inference  of  lung  volumes  from  magnetometer-measured  dimensional  changes 
for  the  general  populus  appears  unlikely  to  be  found. 

Conclusions 

In  the  introduction  we  posed  five  questions  and  feel  we  can  now  comment 
on  them: 

1)  What  dimensional  changes  should  be  measured.  Positions  and  M3 
show  the  most  promise. 

2)  How  are  the  dimensional  changes  related  to  lung  volume  changes. 

The  relationship  is  complex  and  nonlinear,  and  also  exhibits  phase 
shifts  and  hysteresis. 

3-5)  How  reliable  can  the  dimensional  changes  be  measured.  Even  for  a 
given  subject  there  are  breath  to  breath  differences.  Between 
subjects  the  pattern  of  dimensional  changes  is  not  predictable. 

Our  conclusion  is  that  magnetometers  are  useful  for  quantitative 
measurement  of  dimensional  changes  but  because  of  the  complexity  of 
respiration  will  likely  only  have  qualitative  value  in  inferring  lung 
vo  1  ume . 
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Table  1.  Physical  Characteristics  of  Subjects 


SUBJECT 

SEX 

HEIGHT  (cm) 

WEIGHT  (kg) 

CHEST  (cm) 

WAIST  (cm) 

1 

F 

(162) 

(46.3) 

(81.3) 

(58.4) 

2 

M 

(183) 

(72.6) 

(91.4) 

(80.0) 

3 

M 

(190) 

(81.6) 

(96.5) 

(83.8) 

4 

M 

(190) 

(74.8) 

(91.4) 

(78.7) 

5 

M 

(185) 

(76.2) 

(91.4) 

(78.7) 

6 

F 

(162) 

(49.9) 

(86.4) 

(63.5) 

7 

F 

(168) 

(54.4) 

(86.4) 

(66.0) 

8 

M 

(185) 

(74.8) 

(91.4) 

(83.8) 

9 

M 

(183) 

(90.7) 

(114.3) 

(86.4) 

10 

M 

(178) 

(63.5) 

(91.4) 

(78.7) 

11 

M 

(188) 

(90.7) 

(100.3) 

(86.4) 

« 
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TABLE  4  LATERAL  QUIET  BREATHING 


s 

& 

Ml 

Mfi 

s 

1+0 

m5 

.6062+. 51 39 

1+0 

symmetric 

*6 

. 1610+.5795 

.1176+. 5511 

1+0 

m7 

.5803+. 3959 

.4176+. 4808 

.0858+. 5153 

i+0 

Mg 

.4513+. 5254 

.3626+. 4915 

.2598+. 6074 

.4386+. 4985 

1+0 

TABLE  5  LATERAL 

FORCED  BREATHING 

S 

Mi 

M&  Ml 

Mfl 

S 

i+0 

m5 

.5619+.5082 

i+0 

symmetric 

m6 

.3296+. 5982 

.0831+. 6548 

i+0 

m7 

. 3555+. 6198 

.0794+. 6002 

.3194+. 4766  1+0 

Mg 

.5316+. 5390 

.2533+. 5807 

.4172+.4974  .4771+.4415 

1+0 
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Table  6  AP  cummulative  breath  correlation  coefficients 
for  subject  #10.  Note  that  the  values  approach 
the  total  as  the  number  of  breaths  increases. 
Exact  correspondence  is  not  achieved  due  to 
incomplete  breath  portions  at  each  end. 


Number  of 

Breaths _ Ml _ M2 _ M3 _ M4 

1  .9786  .6950  .9636  .6975 

2  .9872  .9542  .9858  .8974 

3  .9691  .8145  .9709  .6934 

4  .9743  .6708  .9699  .6672 

5  .9784  .6184  .9548  .6746 

6  .9802  .5359  .9452  .6918 

7  .9801  .5958  .9472  .7446 

Last  .9799  .6073  .9453  .8184 

TOTAL 

(from  Appendix  2)  .9800  .6327  .9501  .8309 


Table  7 


AP  correlation  coefficients  between  spirometer  and 
magnetometers  for  six  individual  breaths  for  each 
subject.  Notice  significant  variations  only  in  M2 
position. 


SUBJECT  1 


Breath  # 

1 

2 

3 

Mi 

.9943 

.9924 

.9965 

m2 

.7330 

.9153 

.8612 

M3 

.7938 

.7128 

.9286 

m4 

.6046 

.6074 

.6280 

SUBJECT  #2 

Mi 

.9917 

.9944 

.9913 

m2 

-.2067 

-.0136 

-.2203 

M3 

.8693 

.9658 

.8277 

M4 

.9706 

.9940 

.9247 

SUBJECT  #3 

Mi 

.9849 

.9843 

.9921 

m2 

.5853 

.6324 

.8366 

M3 

.9468 

.9757 

.9898 

m4 

.7945 

.9773 

.9880 

SUBJECT  #4 

Mi 

.9975 

.9973 

.9959 

m2 

.8896 

.8783 

.8564 

m3 

.9631 

.9492 

.9458 

m4 

.8987 

.9531 

.9871 

SUBJECT  #5 

Mi 

.9442 

.9234 

.9551 

m2 

.9446 

.9805 

.9271 

m3 

.9562 

.9816 

.9529 

m4 

.9713 

.9590 

.9711 

SUBJECT  #6 

Mi 

.9906 

.9959 

.9953 

m2 

.5938 

-.0588 

.0485 

m3 

.9648 

.7155 

.9594 

m4 

.9629 

.6785 

.9516 

4 

5 

6 

AVG  ♦  SD. 

.9859 

.9869 

.9677 

.9873  +  .010 

.9194 

.9488 

.9455 

.8872  +  .082 

.8702 

.6305 

.8504 

.7977  .109 

.6404 

.8385 

.8791 

.7147  .116 

.9979 

.9876 

.9905 

.9922 

+  .004 

.4364 

-.1368 

.1806 

.0081 

+  .259 

.9765 

.9374 

.9386 

.9192 

+  .058 

.9941 

.9764 

.9812 

.9735 

+  .026 

.9933 

.9948 

.9908 

.9900 

+  .004 

.6728 

.8241 

.9572 

.7515 

+  .143 

.9841 

.9912 

.9937 

.9802 

+  .018 

.9891 

.9918 

.9961 

.9561 

~  .079 

.9933 

•  •  • 

•  •  • 

.9960 

+  .002 

.9247 

•  •  • 

•  •  • 

.8873 

♦  .029 

.9788 

•  #  # 

000 

.9592 

+  .015 

.9572 

•  •  a 

000 

.9490 

+  .037 

.9581 

.9725 

.8983 

.9419 

+  .027 

.9659 

.9546 

.9651 

.9563 

+  .019 

.9908 

.9777 

.9771 

.9727 

+  .015 

.9895 

.9856 

.9897 

.9777 

+  .012 

.9974 

.9972 

.9944 

.9951 

+  .002 

.1482 

.2152 

-.7758 

-.0208 

+  .453 

.9676 

.8804 

.8973 

.8981 

♦  .097 

.9552 

.8816 

.7841 

.8689 

+  .115 
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SUBJECT  #7 


Ml 

.9964 

.9902 

.9885 

.9963 

.9972 

.9982 

.9945  +  .004 

M2 

-.9331 

-.6511 

.2943 

.9136 

.2874 

-.5368 

-.1043  ♦  .710 

m3 

.9697 

.8868 

.8526 

.9957 

.9591 

.9217 

.9309  +  .054 

M4 

.9896 

.9187 

.6794 

.9729 

.9658 

.9857 

.9187  +  .119 

SUBJECT  #8 

Mi 

.9146 

.9708 

.9888 

.9262 

.9693 

.9913 

.9602  +  .032 

m2 

.5116 

.6896 

.7786 

.8107 

.9004 

.9275 

.7697  7  .153 

M-i 

.8628 

.6582 

.8355 

.8398 

.9545 

.9841 

.8558  7  .115 

M4 

.5478 

.4092 

.5868 

.5822 

.9237 

.9571 

.6678  +  .221 

SUBJECT  #9 

Mi 

.9534 

.9224 

.8822 

.9159 

.9576 

.9569 

.9314  +  .030 

m2 

.8300 

.8505 

.9074 

.9395 

.9759 

.9556 

.9099  +  .059 

M-i 

.9707 

.9551 

.9635 

.9754 

.9889 

.9861 

.9733  +  .013 

m4 

.9685 

.9846 

.9899 

.9886 

.9932 

.9715 

.9827  +  .010 

SUBJECT  #10 

Mj 

.9786 

.9863 

.9534 

.9964 

.9905 

.9894 

.9824  +  .015 

m2 

.6949 

.9518 

.8867 

.8838 

.8882 

.8787 

.8640  +  .087 

m3 

.9636 

.9844 

.9856 

.9890 

.9801 

.9732 

.9793  +  .009 

M4 

.6975 

.8892 

.6231 

.9794 

.9729 

.9845 

.8558  ♦  .158 

SUBJECT  #11 

Mi 

.9609 

.9836 

.8344 

.9497 

.9764 

.9377 

.9404  +  .055 

M2 

.9131 

.9914 

.8977 

.9219 

.9497 

.9079 

.9303  +  .035 

«3 

.9456 

.9211 

.9454 

.9701 

.9384 

.9456 

.94*4  +  .016 

«4 

.8217 

.9741 

.8562 

.9602 

.6584 

.9878 

.8704  ♦  .126 
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Table  8  AP  Correlation  coefficients  for  a  single 

forced  breath  for  each  subject.  Notice  no 
remarkable  differences  from  those  of  quiet 
breathing. 


SUBJECT 

Mi 

m2 

m3 

M4 

1 

.9975 

.9782 

.8801 

.9275 

2 

.9829 

.6253 

.9213 

.9537 

3 

.9888 

.4612 

.9698 

.8537 

4 

.9964 

.9564 

.9381 

.8221 

5 

.9706 

.9789 

.9169 

.9198 

6 

.9935 

.2457 

.7463 

.4557 

7 

.9826 

-.7417 

.8894 

.8696 

8 

.9948 

.4384 

.8431 

.9616 

9 

no  measurements  made 

10 

.9844 

.6871 

.8938 

.8866 

11 

.9555 

.7147 

.9204 

.9243 

Average 

.9847+. 01 3 

.5343+. 512 

.8919+. 062 

.8574+. 148 
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Table  9 

AP  Correlation  coefficients  for 
showing  individual  correlations 
(I),  expiration  (E)  and  the  total 

a  single  breath 
for  inspiration 
breath  (T) 

SUBJECT 

TEST 

Mi 

M2 

Jh _ 

"A 

1 

I 

.9964 

.9524 

.9964 

.8491 

E 

.9971 

.9857 

.9919 

.9377 

T 

.9948 

.8823 

.9837 

.8279 

2 

I 

.9989 

.2029 

.9981 

.9978 

E 

.9909 

-.3788 

.9001 

.9784 

T 

.9874 

-.1569 

.9352 

.9756 

3 

I 

.9951 

.9739 

.9945 

.9957 

E 

.9977 

.7081 

.9913 

.9976 

T 

.9948 

.8240 

.9912 

.9917 

4 

I 

.9975 

.9608 

.9853 

.9764 

E 

.9991 

.9433 

.9892 

.9942 

T 

.9973 

.8783 

.9491 

.9531 

5 

I 

.9562 

.9859 

.9944 

.9961 

E 

.9773 

.9672 

.9879 

.9966 

T 

.9551 

.9270 

.9529 

.9711 

6 

I 

.9997 

.8623 

.9989 

.9949 

E 

.9965 

-.4624 

.7004 

.6487 

T 

.9959 

-.0588 

.7155 

.6785 

7 

I 

.9973 

.9598 

.9954 

.9949 

E 

.9978 

.9013 

.9969 

.9869 

T 

.9963 

.9136 

.9957 

.9728 

24 


V 


8 

i 

.9962 

.9822 

.9967 

.9966 

E 

.9965 

.9733 

.9883 

.9869 

T 

.9913 

.9275 

.9841 

.9570 

9 

I 

.9811 

.9311 

.9784 

.9976 

E 

.9982 

.9934 

.9982 

.9971 

T 

.9582 

.9568 

.9865 

.9722 

10 

I 

.9974 

.9479 

.9956 

.9851 

E 

.9981 

.7965 

.9821 

.9845 

T 

.9964 

.8838 

.9890 

.9794 

11 

I 

.9781 

.9699 

.9601 

.9851 

E 

.9957 

.9911 

.9959 

.9960 

8 

T 

.9404 

.9267 

.9544 

.9879 

Table  10  Correlation  coefficients  for  the  single  breath 
for  which  the  fourier  results  (Appendix  3)  are 
given 


SUBJECT 

Mi 

m2 

m3 

m4 

1 

.9965 

.8612 

.9286 

.6280 

2 

.9979 

.4364 

.9765 

.9942 

3 

.9933 

.6729 

.9841 

.9891 

4 

.9975 

.8896 

.9631 

.8987 

5 

.8983 

.9650 

.9771 

.9897 

6 

.9974 

-.1482 

.9675 

.9552 

7 

.9963 

.9136 

.9957 

.9729 

8 

.9913 

.9275 

.9841 

.9571 

9 

.9224 

.8505 

.9551 

.9846 

10 

.9534 

.8868 

.9856 

.6231 

11 


8344 


8977 


9454 


8562 


Figure  3.  A  typical  time  trace  of  the  five  signals  (spirometer, 

M. ,  M_,  M, »  M  )  for  subject  #5  during  a  forced  breathing 
manuever  Is  shown.  Notice  the  maximal  breath  in  the 
eighth  breath  see  text  for  details. 
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Figure 


4.  Hysteresis  curves  of  a  given  breath  for  a  typical  subject. 
Notice  the  more  pronounced  hysteresis  in  the  M-  (bottom) 
position  with  respect  to  the  M,  (top)  position.  (Next 
page  M2(top),  M^(bottom)). 
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Figure  5.  Hysteresis  curves  for  a  given  subject  showing  two  successive 
breaths  at  the  position.  Note  the  variations.  Greater 

variations  were  seen  in  each  of  the  other  three  AP  positions 
on  these  breaths. 
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Figure  6.  Lateral  (M,)  position  hysteresis  curves  showing  significant 
variations  on  successive  breaths  (Three  curves;  continued 
next  page). 
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Appendix  1 


All  data  analysis  and  collection  was  performed  on  a  Mine  11  with  the 
controlling  software  written  in  the  high  level  language,  BASIC.  The  pro¬ 
gram  listings  as  well  as  a  brief  program  explanation  are  given  below.  The 
programs  were  written  so  as  to  "prompt"  the  user  for  the  specific  input  as 
well  as  to  instruct  the  user;  this  to  allow  other  users  to  easily  implement 
the  programs  on  continued  investigation. 

A  few  points  should  be  highlighted  at  this  time.  The  test  data  files 
were  stored  as  one  dimensional  integer  arrays  of  three  thousand  data 
entries.  They  were  stored  as  a  sequence  of  spirometer,  M3,  M2,  M3,  M4  (or 
M5  ,  Mg  ,  M7,  Mg)  repeated  600  times.  Approximately  twenty  four  data  files 
each  of  3000  digitized  points  could  be  stored  per  secondary  (data)  diskette 
(Scotch,  8  inch,  soft  sector,  double  density,  RX02  format). 

In  programs  where  individual  breaths  were  analyzed  (eg.  ANALY  3  .  BAS  ) 
the  program  asks  for  the  spirometer  file  place  (0  -  600)  to  define  the  given 
breath.  This  is  exemplified  in  figure  A.  1 .  Thus  by  inputing  any  max/min 
value  for  the  filplace,  any  breath,  portion  of  a  breath  or  sequence  can  be 
analyzed . 

As  a  final  note,  all  output  data  file  names  for  output  on  disk  drive  2 
must  be  inputed  before  program  execution.  And,  in  the  program  listings, 
1/2  stands  for  and  1/4  stands  for  ">". 

The  programs  are: 

(1)  LUNG  -  The  program  LUNG  controlled  the  data  collection  process  (eq. 
sampling  rate,  sampling  time).  1.  five  simultaneous  signals  were 
converted  to  digitized  form,  stored,  and  graphed  on  the  screen  for  a 
qualitative  "look"  at  the  data. 
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(2)  ANALYP  -  The  program  ANALYP  recalled  the  data  file  from  storage  (from 
LUNG)  and  calculated  the  correlation  coefficients  as  given  in 
Appendix  2.  The  results  were  both  stored  and  displayed  on  the 
screen. 

(3)  RESLT2  -  The  program  RESLT2  recalled  the  data  stored  by  ANALYP  for  any 
number  of  subjects  and  averaged  and  gave  standard  deviations  for  the 
results  (as  given  in  tables  2-5). 

(4)  MCAL  -  The  program  MCAL  allowed  us  to  perform  calibration  tests  on 
each  magnetometer  based  on  separation  distances,  rotational  effects 
and  adjustment  of  relative  gain. 

(5)  SCAL  -  The  program  SCAL  was  used  to  obtain  calibration  curves  for  the 
spirometer  as  output  vs  volume. 

(6)  MAX  4  -  The  program  MAX  4  was  used  to  scan  each  data  file  to  determine 
the  relative  maxima  and  minima  fileplaces  for  use  as  described  above. 
Due  to  the  experimental  nature  of  the  data,  abolute  extrema  were  rare, 
and  results  were  often  doubled  checked  visually  using  MANGE2. 

(7)  MANGE2  -  The  program  MANGE2  was  used  to  visually  inspect  data  and  to 

generate  cross  plots  (hysteresis  curves). 

(8)  ANALY3  -  ANALY3  is  essentially  the  same  as  ANALYP  except  that  one  could 

keyboard  control  which  breath  or  parts  thereof  that  correlations  were 
to  be  calculated.  The  results  of  MAX  4  were  used  here. 

(9)  The  following  "chained  programs"  were  used  to  calculate  the  fourier 

coefficients  (up  to  seven  harmonics)  for  the  given  breath.  They 
allowed  visual  inspection  of  the  breath,  calculation  of  the  coeffi¬ 
cients,  visual  comparison  of  the  results  (see  Figure  8),  normaliza¬ 
tion  of  results  with  respect  to  the  spirometer  and  storage. 
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(i)  FMAIN 


(ii) 

FGRAPH 

(iii) 

FCOEFF 

(iv) 

FCOMP 

(v) 

FNORM 

(vi) 

FOUT 

(vii ) 

NS  TO  RE 

LUNG.BAS 


10  BJ«'  ****  ******  ******  ******  ******  ******' 

20  DIM  M%(3000)  T«0 

30  DISPLAY  CLEAR 

40  PRINT  ■•*•*••**•*******••******•••*****•*****• 

41  PRINT  '*  This  program  collects  data  from  *' 

42  PRINT  '*  channels  0-4  and  writes  to  a  file  *' 

43  PRINT  ’*  or  to  the  screen.  •' 

44  PRINT  •**********e*e*e*e**e****e*ee**e*eeee*e' 

45  PRINT  " 

50  PRINT  \  PRINT  'Enter  the  filename  for  output  ';  S  UN  PUT  F$ 

55  PRINT  'Enter  the  file  number  1-10';  \  INPUT  F5 

56  PRINT  'Enter  the  number  of  data  points  to  be  taken';  \  INPUT  D1 

58  PRINT  "  PRINT  " 

59  PRINT  'notes  did  you  substitute  filename  in  #200' 

60  PRINT  \  PRINT  \  PRINT  'Are  you  ready  KRETURNH’;  \  INPUT  A* 

70  PRINT  \  PRINT  'input  to  begin  after  beep.....' 

80  FOR  J»1  TO  100  \  NEXT  J 
90  1*1 

100  SCHEDULECINTERVAL',  1,140) 

110  PRINT  CHR$(7) 

120  AIN(,M%(),D1,.075,0,5) 

130  IF  1=0  THEN  150  \  GO  TO  120 
140  1=0  RETURN 

150  PRINT  CHR$(7)  \  PRINT  \  PRINT  Input  from  channel  0-4  now  complete.' 

160  GRAPH(,„M%()) 

165  PRINT  '• 

170  PRINT  'Output  to  disk  or  to  screen  or  do  nothing  (D  or  S  or  N)';  \  INPUT  A$ 

190  IF  AS='D'  GO  TO  200  \  IF  AJ-'S'  GO  TO  230  \  IF  AS-'N'  GO  TO  270  \  GO  TO  170 

200  OPEN  'DKlsF.DAT'  FOR  OUTPUT  AS  FILE  #F5 

210  FOR  K=0  TO  Ol-l  \  PRINT  #F5,M%(K)  \  NEXT  K 

220  CLOSE  #F5  \  DISPLAY  CLEAR  \  GO  TO  270 

230  DISPLA Y_C LEAR 

240  FOR  K=0  TO  Dl/S-1  \  K%=5*K 

250  PRINT  USING  BJ,K%,M%(X%),M%(K%+l),M%(K%+2),M%(K%+3),M%(K%+4) 

260  NEXT  K 
270  PRINT 

290  DISPLAY  CLEAR \  STOP  \  END 
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ANALYP.BAS 


10  DIM  M%(4,599),S(4),Q(4),T(414),Cl4,4) 

20  R*0  ZJa'.DAT  WJ*'#################' 

28  PRINT  '•••••**e*********e****ee**ee**see- 

29  PRINT  '*  This  program  calculates  the  *' 

30  PRINT  '*  sum,  the  sum-squared  and  the  *' 

31  PRINT  '*  correlation  coefficients  for  •' 

32  PRINT  '*  channels  0-4  *' 

33  PRINT  '*•**••***•**••*******•**•*•••**«• 

34  PRINT  "  PRINT  " 

40  PRINT  PRINT  'Enter  the  filename  for  input  \  INPUT  F$ 

45  PRINT  'Enter  the  file  number  1-10';  INPUT  F6 

4b  PRINT  'Enter  the  number  of  data  points  collected';  \  INPUT  Cl 

47  PRINT  PRINT 

48  PRINT  '  note:  did  you  substitute  filename  into  #60' 

49  C 1*3000 

50  OPEN  FJVO.DAT  FOR  OUTPUT  AS  FILE  #10 
60  OPEN  'DK1-.F.DAT  FOR  INPUT  AS  FILE  #F6 

70  FOR  J=0  TO  Cl/5-1  \  FOR  1*0  TO  4  \  INPUT  #F6,M%(I,J)  \  NEXT  I  \  NEXT  J 
80  CLOSE  #F6 
90  FOR  1*0  TO  4  \  Sd)=0 
95  L=C  1/5-1 

100  FOR  K*0  TO  4  \  T(I,K)*0  \C(I,K)*0  \  NEXT  K 
110  FOR  J*0  TO  L  \  Sa)»S(D+M%a,J)  \  NEXT  J 
115  FOR  K=0  TO  I 

120  FOR  J=0  TO  L  \  PRINT  \  Td,K>*T(I,K)+l*M%a,J)*M%(K,J)  \NEXT  J 
125  NEXT  K 

130  FOR  K*0  TO  I  \  Ca,K>-Ta,K)-S(I)*S(K)/(L+l)  \  NEXT  K  \  NEXT  I 
135  FOR  1*0  TO  4 

140  Qd)*SQR(ABS(ca4)))  \  FOR  K*0  TO  I  \  Ca,K)*C(I,K)/QCI)  \  NEXT  K 

145  NEXT  I 

146  FOR  1*0  TO  4 

150  FOR  K-0  TO  I  \  CO,K)*INT(Ca,K)/Q(K)*lOOOO+.5)/10OOO  \  NEXT  K 

151  NEXT  I 

160  FOR  1*0  TO  4  \  FOR  K-I  TO  4  \  Td,K)*T(K®  \  Cd,K)=C(K4)  \  NEXT  K  \  NEXT 
165  PRINT  \  PRINT  WJ  \  PRINT  *#’ 

170  PRINT  '♦  'StF$J*STR$(R)J*Z$  \  PRINT  '#'  \  PRINT  WJ 
190  PRINT  \  PRINT  '***SUM***'  \  PRINT  S(0),S(1),S(2),S(3),S(4) 

205  PRINT  \  PRINT  '***INNER  PRODUCT***' 

210  FOR  1*0  TO  4  \  PRINT  Ta,0),TCI,l),Ta,2),Ta,3),TCI,4)\  NEXT  I 
225  PRINT  \  PRINT  ’***CORRELAT10N  COEFFICIENT***’ 

230  FOR  1-0  TO  4  \  PRINT  Cd,0),C(I,l),CCI,2),Ca»3>,C(I,4)  N  NEXT  I 

240  FOR  1*0  TO  4  \  FOR  J*0  TO  4  \  PRINT  #10,CCI,J)  \  NEXT  J  NEXT  I 

245  CLOSE  #10 

246  PRINT  'AGAIN';  \  LINPUT  H$ 

247  IF  HS»'A'  THEN  PRINT  '***' 

250  DISPLAY  CLEAR  \  STOP  \  END 
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RESLTZ.BAS 


10  DIM  S(4,4),C(4,4),A(4,4),Sl(4t4),D(4,4) 

15  FOR  K«0  TO  4  \  FOR  L-0  TO  4  \  S(K,L)-0  \  NEXT  L  \  NEXT  K 

17  N-0 

18  PRINT  'Enter  the  output  filename';  \  LINPUT  G$ 

Z0  PRINT  'Enter  the  input  filename';  \  LINPUT  F$ 

Z5  N-N+l 

30  OPEN  FSA'O.DAT*  FOR  INPUT  AS  FILE  #5 
40  FOR  I«0  TO  4 
45  FOR  J-0  TO  4 
50  INPUT  #5,Cd,J) 

55  NEXT  J 
60  NEXT  I 
70  CLOSE  #5 
80  FOR  1-0  TO  4 
85  FOR  J-0  TO  4 

ioo  sa,j)>sa.J)ca,J) 
io5  sia,J)-sicu)+cn,J)*2 

110  NEXT  J 
1Z0  NEXT  I 

170  PRINT  'Again?';  \  LINPUT  AJ 
180  IF  A8>'Y‘  THEN  20 
ZOO  FOR  K«0  TO  4 
210  FOR  L-0  TO  4 
Z2C  A(1C,L)-S(K,L)/N 

ZZ5  D(KJ.)-SQR(ABS((S  1(K,L)-N*A(K,L)a  2)/(N- 1))) 

230  NEXT  L 
240  NEXT  K 
24*  PRINT  "  PRINT  " 

245  PRINT  'Wish  to  store  this:';  \  LINPUT  Q$ 

246  IF  Q8-'N'  THEN  271 

250  OPEN  GJfc'.QAr  FOR  OUTPUT  AS  FILE  #3 

255  PRINT  *3, 'THE  AVERAGE  VALLES  ARE:' 

256  PRINT  43,'  ' 

260  FOR  1-0  TO  4  \  PRINT  43,Aa,0),ACI,l>,Aa,2),A(U>,ACl,4)  \  NEXT  I 

262  PRINT  #3,'  ' 

263  PRINT  43, THE  STANDARD  DEVIATIONS  ARE: 

264  PRINT  43,'  ’ 

265  FOR  J-0  TO  4  \  PRINT  43,D(JI0),D(J,1),D(J,2),D(J,3),D<J,4)  \  NEXT  J 

270  CLOSE  43 

271  DISPLAY  CLEAR  \  PRINT  AVERAGE  VALUES  ARE:'  \  PRINT 

272  FOR  1-0  TO  4  \  PRINT  A(I,0>,Aa,l),ACU>,A(U>,A!I,4)  '  NEXT  I 

273  PRINT  "  \  PRINT  'THE  STANDARD  DEVIATIONS  ARE:’  \  PRINT  " 

274  FOR  J-0  TO  4  \  PRINT  D(J,0),D<J,l>,DtJ,Z),DW,3>,DU,4)  \  NEXT  J 
280  STOP  \  END 
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MCAL.BAS 


1  REM  DleSEPARATION  DISTANCES,  Sl-AVG.  VOLTAGE/DISTANCE 

2  REM  A3*MEAN  DISTANCE, A4-MEAN  VOLTAGE 

3  REM  B3>SLOPE,  B4-INTERCEPT 

10  REM  THE  PROGRAM  NAME  IS  MCAL 
20  REM  ENTER  CHANNEL  #5  TO  STOP 
SO  DIM  Dl(20), S1(20),V(100),V1(100) 

55  DIM  S(100),A1(20),A2(20),B1(20),B2(20) 

60  1*0  \  M*20 

62  PRINT 

63  PRINT  '*  This  program  allow*  calibration  of  the  ** 

64  PRINT  '*  magnetometer*  eoltage  es.  separation  die-  *' 

65  PRINT  '*  tance...  a  least  squares  fit  can  then  be  *' 

66  PRINT  '•  computed.  •* 

67  PRINT  'eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee** 

68  PRINT  "  \  PRINT  " 

80  FOR  K*0  TO  M\  S1(K)*0  \  D1(K)*0  \  NEXT  K 
86  FOR  K*0  TO  99  \  V(K)-0  \  V1(K)*0  \  NEXT  K 

100  1*1+1  \  J-0  \  Al(l)-0  \  A2(l)*0  \  B1(1)*0  \  B2(l)-0  \  B3*0  \B4«0 

101  PRINT  'ENTER  THE  MAGNETOMETER  NUMBER'*,  \  INPUT  II 

102  I>I1 

103  IF  1*5  THEN  275 
105  DISPLAY  CLEAR 

109  PRINT  •***•*******************•••***' 

110  PRINT  USING  'SET  UP  MAGNETOMETER  PAIR  #  .'4 

111  PRINT  '••***•••**•••*•*•#••••**•••**' 

114  PRINT  '•  PRINT  " 

120  PRINT  'ANOTHER  SEPARATION  DISTANCE  Y  or  N’;  \  UNPUT  Q* 

121  PRINT  " 

130  IF  Q$-*NO'  THEN  232  \  IF  QJ-'N'  THEN  232 

131  J*J+1 

140  PRINT  'PLEASE  GIVE  SEPARATION  DISTANCE  (inches)'!  \  INPUT  D 
150  D1(J)*D 

160  SCHEDULEONTERVAL',  10,200) 

170  AIN(’DISPLAY,,V0,100,.1,I) 

180  IF  L*0  THEN  205 
190  GO  TO  170 
200  L*0  \  RETURN 
nr^PT  av  n  par 

210  FOR  K*0~ TO  99  \  V1(K)«V(K)  \  S(K*1)^(K)+V1(K) 

215  NEXT  K 
220  S1(J)*S(K+1)/K 
225  M*J 
230  GO  TO  120 

232  PRINT  H 

233  PRINT  '  DISTANCE  AVG.  VOLTAGE* 

234  PRINT  " 

235  FOR  L*1  TO  J  \  PRINT  D1(L),S1(L)  \  NEXT  L 

236  PAUSEfi) 
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240  PRINT  "  PRINT  " 

265  PRINT  'DO  YOU  WISH  A  LEAST  SQUARES  ANALYSIS; \  UNPUT  LJ 

270  IF  LJ-'YES'  THEN  GOSUB  300 

271  IF  LS-’Y'  THEN  GOSUB  300 

272  PRINT  *DO  YOU  WISH  TO  STORE  THIS;’!  \  UNPUT  D$ 

273  IF  DJ-'V  THEN  GOSUB  480  \  IF  DS-’YES'  THEN  GOSUB  480 

274  DISPLAY  CLEAR 

275  IF  IMS  THEN  80  \  IF  1-5  THEN  STOP 

300  REM  THIS  IS  THE  LEAST  SQUARES  SUBROUTINE 
310  FOR  K-l  TO  M 
320  A 1  (K+ 1)» A 1  (K)*D  1  (K) 

330  A2(K+1)-A2(K)+S1(K) 

340  NEXT  K 

350  A3*Al(M+l)/M 

360  A4«A2(M+1)/M 

370  FOR  L-l  TO  M 

380  B 1  (L+ 1)  -B 1  (L)+ (D 1  (L)- A3)  *  (S 1  (L)- A4) 

390  B2(L+1)»B2(LWD1(L)-A3)  2 

400  NEXT  L 

410  B3-B1(M+1)/B2(M+1) 

420  B4-A4-B3*A3 
425  PRINT 

430  PRINT  ’THE  SLOPE  IS;';  \  PRINT  B3 
440  PRINT  'THE  INTERCEPT  IS:1;  \  PRINT  B4 
450  PRINT  " 

465  PAUSE(IS) 

466  DISPLAY  CLEAR 
470  RETURN- 

480  REM  THIS  IS  A  SUBROUTINE  TO  STORE  DATA 
490  PRINT  ’ENTER  THE  DATA  FILE  NAME:';  \  UNPUT  SJ 
500  OPEN  SSlr'.DAT  FOR  OUTPUT  AS  FILE  *3 
510  FOR  K-l  TO  J\  PRINT  #3,D1(K),S1<K)  \  NEXT  K 
520  CLOSE  #3 

530  PRINT  \  PRINT  \  PRINT  \  PRINT 
540  RETURN 
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SC  ALBAS 


SO  REM  PROGRAM  NAME  IS  SCAL 

100  DIM  AUOO),Am(100),A2(U),SUOO),S2{ll),Sl(ll),Zl<ll) 

120  PRINT  '•••••••••••••••••••••••••••••••••••••v****************' 

121  PRINT  '•  *' 

122  PRINT  '*  Thu  program  calibrate!  input  voltages  in  terms  *' 

123  PRINT  '*  of  volume  for  the  spirometer.. .begin  at  voL*0  *' 

124  PRINT  '•  *' 

125  PRINT  •*•••••*•••****•***••*****••*******•**•****•***••***•••’ 
130  REM  DESIGNED  TO  INPUT  10  VOLUMES  AND  CALIBRATE  (0-10) 
140  J»0  \  I»1  \  S(0)«0  \  S2(0)*0  \  S1(0)=0  \  Zl(0)-0 

150  SCHEDULE(TNTERVAL’,10,190) 

160  AIN fDISPLA  Y’,A I960, 100,.  1 ,0,) 

180  IF  I»0  GO  TO  190  \GO  TO  160 
190  WO 

200  FOR  K»0  TO  99  \  A(K)*A1%(K)  \  PRINT  A(K)  \  NEXT  K 
210  FOR  L«0  TO  99  \  S(L*1)*S(L)+A(L)  \  NEXT  L 
220  A2(J)*S(100)/100 
225  DISPLAY  CLEAR 

230  PRINT  ’THE  VOLUME  MEASURED  WAS:’;  \  PRINT  J 
240  PRINT  ’AVERAGE  VOLTAGE  WAS:  \  PRINT  A2(J) 

243  PRINT  " 

244  PRINT  ’•••YOU  HAVE  TEN  SECONDS  TO  INCREASE  VOLUME***’ 

245  PAUSE(IO) 

246  PRINT  CHRS(7) 

250  1*1  \  J-J+l  \  S(0)-0 

260  SCHEDULECINTERVAL*. 10,190) 

270  IF  JH*11  THEN  285 
280  RETURN 

285  PRINT  VOLUME  VOLTAGE’ 

286  FOR  I-O  TO  10  \  PRINT  LA  2d)  \  NEXT  I 

287  PRINT  \  PRINT 

290  REM  NOW  HAVE  VOLTAGES  VS  VOLUME  DATA 

300  REM  NEXT  PERFORM  LEAST  SQUARES  DATA  REDUCTION 

310  FOR  J*0  TO  10  \  S2(J+1)=S2(J)+AZ(J)  \  NEXT  J 

320  V3*SZ(11)/U 

330  V4*5 

340  FOR  K*0  TO  10  \  A(K)»(A2(K)-V3)*(K-V4)  \  NEXT  K 
350  FOR  K*0  TO  10  \  S1(K+1)*A(K)+S1(K)  \  NEXT  K 
360  FOR  K-0  TO  10  \  2(K)*(K-V4)/' 2  \  NEXT  K 
370  FOR  K*0  TO  10  \  Zl(K+l)*Z(K)+Zl(K)  \  NEXT  K 
375  K*ll 

380  M*S1(K)/Z1(K) 

390  B*V3-M*V4 

400  REM  VOLTS  *VOLUME*M+B  ie.  UNEAR  RELATIONSHIP 
410  REM  USE  EQN  VOL*(VOLTS-B)/M 

416  PRINT  ’THE  LEAST  SQUARES  INFORMATION  IS:’ 

417  PRINT  \  PRINT 

420  PRINT  THE  SLOPE  IS:’j  \  PRINT  M 
430  PRINT  ’THE  Y  INTERCEPT  IS:’;  \  PRINT  B 
440  PAUSE(20) 

450  DISPLAY  CLEAR 
460  STOP  \  END 


MAX4.BAS 


1  REM 

100  REM  This  program  allows  one  to  {lad  the  relative  extremum 

101  REM  values  in  the  spirometer  data  ret 
120  DIM  M%(3000),S%(600),C1(5Q),C2(SO) 

1Z5  DIM  C3(50),C4(50) 

130  B»0\  D*0 

U0  OPEN  *DK1:  JDH1.DAT*  FOR  INPUT  AS  FILE  #3 
150  FOR  I»0  TO  Z999 
160  INPUT  #3,Mtt0) 

170  NEXT  I 
180  CLOSE  #3 
190  REM 

ZOO  FOR  K»0  TO  599  \  K%«5*K 
Z10  S%(K)»M%(K%) 

ZZ0  NEXT  K 
Z30  REM 

Z40  FOR  K*3  TO  595 

Z60  T1-S9MK)  \  T2»S%(K+1)  \T3*S%(K+2)  \  T4-S%(K-3)  \  T5«S%(K+3) 

270  IF  T1H-T2  THEN  IF  T3H-T2  THEN  IF  T4H-T2  THEN  IF  TSN-T2  THEN  B«B+1  \  C1(B)*K+ 
1 

280  IF  T1H*T2  THEN  IF  T3H-T2  THEN  IF  T4M-T2  THEN  IF  T5H-TZ  THEN  D=D+1  \  CZ(D)»K+ 
1 

290  NEXT  K 

300  IF  BHD  THEN  J1*B 

301  IF  DH-B  THEN  J1*D 

305  PRINT  ’  FILE  PLACE  MAX  FILE  PLACE  MIN' 

306  PRINT  \  PRINT 

310  FOR  J-l  TO  J1  \  PRINT  C1(J),S%(C1(J)!,C2(J),S%(C2(J))  \  NEXT  J 
320  REM 

325  PRINT  PRINT 

330  PRINT  ’WISH  TO  STORE  MAX  /  MIN  VALUES';  \  LINPUT  S$ 

340  IF  S$»’N'  THEN  380 

345  PRINT  PRINT  'ENTER  THE  FILENAME  PLEASE';  \  LINPUT  F$ 

350  OPEN  FSfc'.DAT  FOR  OUTPUT  AS  FILE  #4 

355  PRINT  *4,'  FILE  PLACE  MAX  FILE  PLACE  MIN' 

360  PRINT  *4, '•*****•*•**•••»•••**•**•••*•*•••*••*•*•*•**•**•' 

365  FOR  J»1  TO  J1  PRINT  ♦4,C1(J),S%(CKJ)),C2(J),S%(C2(J))  \  NEXT  J 

370  CLOSE  #4 

380  PRINT  \  PRINT 

390  STOP  \  END 
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MANGE2.BAS 


I  REM  Program  name  HELLO  or  MANAGE 

5  REM  THIS  PROGRAM  ALLOWS  ONE  TO  OPEN  A  MAGNETOMETER  DATA 

6  REM  FILE  AND  ...  GRAPH, TRANSFER  OR  DISPLAY  THE  DATA 

7  REM  note:  must  tub  filename  into  #30  and  #120 

8  Bt>'  ####  #####  #####  #####  #####  #####' 

10  DIM  M9M3000) 

II  DIM  S%(600),  Y%(600) 

30  OPEN  'DKlsF2.DAT  FOR  INPUT  AS  FILE  #3 
40  FOR  I=*0  TO  2999 
50  INPUT  #3,M%® 

60  NEXT  I 
70  CLOSE  #3 

90  PRINT  'DO  YOU  WISH  TO  COPY  FILE  TO  ANOTHER  DISK';  \  LINPUT  F$ 

100  IF  F$*'N'  THEN  163 
105  DISPLAY  CLEAR 

110  PRINT  *PLACE  NEW  DISKETTE  INTO  DRIVE  2' 

115  PAUSEU5)  \  PRINT  \  PRINT 

120  OPEN  'DKlsF.DAT*  FOR  OUTPUT  AS  FILE  #5 

130  FOR  1*0  TO  2999 

140  PRINT  #5,M%0) 

150  NEXT  I 
160  CLOSE  #5 

163  PRINT  \  PRINT  \  PRINT 

196  PRINT  'GRAPH  THE  RESULTS';  \  LINPUT  GS 

198  IF  G S=W  THEN  205 

199  GRAPH(„,M%0) 

205  DISPLAY  CLEAR  \  PRINT  'GRAPH  SPIROMETRIC  RESULTS';  \  LINPUT  VS 

210  IF  VJ-'Y'  THEN  GOSUB  325 

211  PRINT  'DISPLAY  THE  NUMERICAL  RESULTS’;  \  LINPUT  D$ 

212  IF  DS-’Y'  THEN  GOSUB  385 

215  PRINT  'WISH  TO  GRAPH  HYSTERESIS’;  \  INPUT  HS 

216  IF  HS*'N'  THEN  320 

217  PRINT  'ENTER  DESIRED  MAGNETOMETER  FOR  COMPARISION’;  \  INPUT  I 

218  PRINT  'ENTER  THE  FILEPLACES  DEFINING  DATA  RANGE*;  \  INPUT  M1,M2 
230  FOR  K*M1  TO  M2  \  K%«5*K 

240  S%(K)»M%(K%)  \  Y%( K)  *M%( K%+I) 

250  NEXT  K 

260  GRAPH(„Y%(),S%0) 

300  PRINT  'AGAIN';  \  UNPUT  AS 
305  DISPLAY_CLEAR 
310  IF  AS»'Y'  THEN  196 

315  PRINT  'WISH  TO  STORE  THIS';  \  UNPUT  SS 

316  IF  SS-'Y'  THEN  GOSUB  500 
320  STOP 

325  REM  SUBROUTINE 

326  PRINT  'ENTER  CHANNEL  NO.  0-4';  \  INPUT  M7 
330  FOR  J5-0  TO  599 

335  S%(J5)«M%(J5*5+M7) 

340  NEXT  J5 
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3S0  GRAPH(,„S%0) 

360  LABEL(, 'SPIRO METRIC  CURVE* ,V) 

370  PAUSE(IS)  \  OISPLA Y_C LEAR 

380  RETURN 

38S  REM  SUBROUTINE 

390  FOR  K«0  TO  S99  \  KW»5*K 

400  PRINT  USING  B$,K,M%<K%),M%(K%+1),M%(K%+  2),M%(K%+3),M%(K%>4) 

410  NEXT  K 

420  DISPLAY_CLEAR 

430  RETURN 

500  PRINT  ’ENTER  OUTPUT  FILENAME*;  \  UNPUT  F* 

510  OPEN  FSfc'.DAT  FOR  OUTPUT  AS  FILE  #3 
520  FOR  I*M1  TO  M2 
530  PRINT  #34,S%(I),Y%d) 

540  NEXT  I 
550  CLOSE  #3 
560  RETURN 


ANALY3.BAS 


100  REM  Thi*  program  calculate*  correlation  coefficient* 

110  REM  for  inspiration/ expiration  parts  of  the  data 
ZOO  DIM  M%(3000),T*(525),E%(525) 

210  DIM  M1(Z5),C(2),MZ(2S) 

220  REM 

230  Z$*'CHANNEL  INSPIRATION  EXPIRATION* 

240  OPEN  *DKl:F.DAr  FOR  INPUT  AS  FILE  *3 

250  FOR  1*0  TO  2999 

260  INPUT  #3,M%(D 

270  NEXT  I 

280  CLOSE  #3 

290  PRINT  'Enter  the  number  of  MAX, MIN  value*';  \  INPUT  B 
300  FOR  J*1  TO  B 

310  PRINT  ^nter  the  MAX, MIN  fUeplace';  \  INPUT  M1(J),M2(J) 

320  NEXT  J 
330  A 1*0 

340  FOR  K*1  TO  B-l 

350  FOR  C*M2(K)  TO  M1(K+1) 

360  I%(A1)*M%(C*5) 

370  A1*A1+1 
380  NEXT  C 
390  NEXT  K 
400  A1*A1-1 
410  A2*0 

420  FOR  K-l  TO  B-l 
430  FOR  L-M10Q  TO  M2(K) 

440  E9t<A2)*M%<L*5) 

450  AZ-A2+1 
460  NEXT  L 
470  NEXT  K 
480  A2-A2-1 
490  DISPLAY  CLEAR 

500  PRINT  ’WISH  TO  SEE  EXPIRATION/INSPIRATION  GRAPHS*;  \  LINPUT  El 
510  IF  E8»’N*  THEN  560 

520  REGK>NCUPPER*,1)  \  REGION  CLOWER*. 2) 

530  GRAPH(,„I%<),2„1)  \  GRAPH(,„E%(),2„2) 

540  LABEL(,'EXPIRATTON','.',2)  \  LABEL  (.INSPIRATION','.',  1) 

550  PAUSEG5)  s  DISPLAY  CLEAR 
560  PRINT  \  PRINT 

570  PRINT  'WISH  TO  CALCULATE  CORRELATION  COEFFICIENTS';  \  LINPUT  C$ 

580  IF  C$*'N'  THEN  1020 

590  DISPLAY  CLEAR 

600  PRINT  Z1 

610  PRINT  \  PRINT 

620  FOR  M7*0  TO  4 

630  Pl-0  \  S1»0  \  S2*0  \  S3*0  \  S4-0  \A3-0  \  A 4*0  \  D1*0  \  D2*0  \  N-0 

640  C(l)-0  \  C(2)*0 

650  FOR  K»1  TO  B-l 

660  FOR  C*M2(X}  TO  M1(K+1) 

670  N-N+l 
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640  X»5*C  \  Y»5*OM7 
690  S1-S1+M%(X) 

700  S2*S2+M56(Y) 

710  S3*S3+M%(X)«2 
720  S4«S4+M%<YW 
730  Q1=MH(X)\  Q2»M%(Y) 

740  P1»P1+Q1*Q2 
750  NEXT  C 
760  NEXT  K 
770  A3-S1/N 
780  A4-S2/N 

790  D1«SQR(ABS(1/(N-1)*(S3-N*A3a2))) 

800  D2=SQR(ABS(1/(N-1)*(S4-N*A4a2))) 

810  C(1)»(1/(N-1))*((P1-N*A3*A4)/(D1*D2)) 

820  P2»0  \  R1»0  \  R2»0  \  R3*0  \  R4-0  \  A5»0  \  A6»0  \  D3»0  \  D4»0  \  M=0 

830  FOR  K»1  TO  B-l 

840  FOR  OMl(K)  TO  M2(K) 

850  M-M+l 

860  X»5*C  \  Y=5»C+M7 
870  R1»R1+M?WX) 

880  R2»R2+M94(Y) 

890  R3«R3+M%(X)»2 
900  R4»R4+M%(Y)*  2 
910  Q3-M94(X)  \  Q4-M%(Y) 

920  P2*P2*Q3*Q4 
930  NEXT  C 
940  NEXT  K 
950  A5-R1/M 
960  A6-R2/M 

970  D3«SQR(ABS(1/(M-1)*(R3-M*A5«2)» 

980  D4“SQR(ABS(1/(M-1)*(R4-M*A6A  2))) 

990  C(2W1/(M-1))*((P2-M*A5»A6)/(D3*04)) 

1000  PRINT  M7,C(1),C(2) 

1010  NEXT  M7 
1020  STOP  \  END 
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F MAIN  .BAS 


100  REM  Program  umt  la  FMAIN 
120  COMMON  M%(3000) 

140  OPEN  ’DKliF2.DAT  FOR  INPUT  AS  FILE  *3 
160  FOR  I«0  TO  2999 
180  INPUT  »3,M9ta) 

200  NEXT  1 
220  CLOSE  #3 

240  PRINT  ’Wi»h  to  »e«  spirometric  curve’;  V  UNPUT  S3 
260  IF  SS-'Y'  THEN  CHAIN  'FGRAPH.BAS' 

280  PRINT  'Wiah  to  calculate  fourier  coefficient*';  \  LINPUT  F$ 
300  IF  F$»'Y'  THEN  CHAIN  TCOEFF.BAS' 

320  PRINT  'Wish  to  print  out  FOURIER  RESULTS?’;  \  LINPUT  R3 
340  IF  RS-’Y’  THEN  CHAIN  'FOUTJAS' 

360  STOP  X  END 


« 

FGRAPH.BAS 

100  REM  Program  FGRAPH.BAS  to  be  used  with  FMAIN.BAS  and  FCOEFF  in  chain 
120  DIM  S%(600) 

140  COMMON  M 94(3 000) 

160  PRINT  'Which  channel  please';  \  INPUT  C 
170  FOR  J»0  TO  599  \  S%(J)»0  \  NEXT  1 
180  FOR  1=0  TO  599 
200  S%0)=M%(I*5+C) 

220  NEXT  I 

240  GRAPH(,„S%()) 

260  PAUSE(S)  \  DISPLAY  CLEAR 
280  PRINT  'AGAIN';  \  LINPUT  AS 
300  IF  AS-'Y'  THEN  160 
320  CHAIN  'FCOEFF .BAS' 
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FCOEFF.BAS 


100  REM  program  FCOEFF.BAS 
120  DIM  A(5,7),B(5,7) 

140  COMMON  M%(3000),Al(4)1G(5,7)fP(5,7),N,N2,Ml,M2 

160  PRINT  'Enter  number  of  harmonica  to  be  found';  \  INPUT  N2 

180  PRINT  'Enter  max  datafile  place*  to  define  data  range';  \  INPUT  Ml, M2 

200  N»M2-M1+1 

220  FOR  K»0  TO  4\  FOR  T»M1  TO  M2  \  A1(K)«A1(K)+M%(T*5+K)  \  NEXT  T  \  NEXT  K 
240  FOR  K*0  TO  4  \  A1(K)*A1(K)/N  \  NEXT  K 
246  FOR  J»1  TO  N2 

250  FOR  K»0  TO  4  \  A(K,J)»0  \  B(K,J)»0  \  NEXT  K 

255  NEXT  J 

260  FOR  J»1  TO  N2 

280  FOR  K»0  TO  4 

300  FOR  T»M1  TO  M2 

320  A(K,J)*A(KfJ)+M%(T*5+K)*COS(2»FT*J»T/N) 

340  B(K,J)»B(K,J)+M%(T*5*K)*SIN(2»PI*J*T/N) 

360  NEXT  T 

380  NEXT  K 

400  NEXT  J 

440  FOR  J-l  TO  N2 

460  FOR  K-0  TO  4 

480  A(K,J)=A(K,J)*Z/N 

500  B<K,J)»B(K,J)*2/N 

520  G(K,J)*SQR(A(K,J)a  2*B(K,J)  a2) 

540  P(K^)«ATN(A(K^X)/B(K,J)) 

560  IF  B(K,J)H0  THEN  P(K,J)-P(K,JHPI 
580  NEXT  K 
600  NEXT  J 

620  CHAIN  TCOMP.BAS' 
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FCOMP.BAS 


100  REM  program  FCOMP.BAS 
110  DIM  S%(600),F%<600) 

1Z0  COMMON  M%(3000),A1(4),G(5,7),P(5,7),N,NZ,M1,MZ 

130  PRINT  'Wish  to  compare  the  results  with  original  file';  \  LIN  PUT  C  J 

135  IF  C$*'N‘  THEN  540 

136  DISPLAY_CLEAR 

140  PRINT  'Enter  the  channel  *  lot  comparison:1;  \  INPUT  M7 
160  FOR  T*M1  TO  MZ  \  F%(T)=0  \  NEXT  T 
180  FOR  T=M1  TO  MZ 
ZOO  FOR  J»1  TO  NZ 

ZZ0  F%(T)=F%(T)+G(M7,J)*SIN(Z*PI*T*J/N+P{M7,J)) 

Z40  NEXT  J 

Z60  F%(T)=F%(T)+A  1  (M7) 

Z80  S%(T)=M%(T*5+M7) 

300  PRINT  S%(T),F%(T) 

3Z0  NEXT  T 

340  DISPLAY  CLEAR 

360  REGION  CUPPER', 1)  \  REGION  ('LOWER',!) 

380  GRAPH(,„S%(),Z„1)  \  GRAPH(,„F%(),2„Z) 

400  LABEL (,'SPIRO M ETRIC  RESULT  one  breath','.’,l) 

4Z0  LABEL(, 'FOURIER  ANALYSIS  REPRESENTATION','.',!) 

440  PRINT  'The  number  of  harmonics  found  is:',NZ 

460  PRINT  \  PRINT  'CLEAR  THE  SCREEN?';  \  LINPUT  QS 

480  PRINT  'AGAIN';  \  LINPUT  AS 

500  IF  AS='Y'  THEN  136 

SZ0  IF  QS='Y*  THEN  DISPLAY  CLEAR 

540  CHAIN  'FOUT.BAS* 
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FNORM. BAS 

100  REM  program  FNORM. BAS 

120  COMMON  M*(3000),A1(4)fG<5,7),P(5,7),N,N2,M1,M2 

140  PRINT  'Wish  to  print  out  normalized  fourier  results';  \  LINPUT  R$ 

180  DISPLAY  CLEAR 
200  IF  Rgs'JT'  THEN  510 

220  PRINT  '•*•  FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA 
240  PRINT  \  PRINT 

260  PRINT  '  NORMALIZED  AMPLITUDES' 

280  PRINT 

300  PRINT  'SPIROMETER  MAGI  MAG2  MAG3  MAG4 ' 

310  G1=G(0,1) 

320  FOR  Jsl  TO  N2 

341  PRINT  G(0,J)/G1 ,G(1 ,J)/G(1 ,1) ,G(2,J)/G(2, 1) ,GC 3 , J)/G(3 , 1 ) ,G(4 , J)/G(4 , 1 ) 

342  REM  PRINT  G(0 , J) /G 1 ,G( 1 , J ) /G 1 ,G< 2 , J ) /G 1 ,G( 3 , J) /G 1 ,G(4  ,  J ) /G 1 
360  NEXT  J 

380  PRINT 

400  PRINT  '  NORMALIZED  PHASE  ANGLES' 

420  PRINT 

440  PRINT  'SPIROMETER  MAGI  MAG2  MAG3  MAG4 ' 

450  PI sPI/2 

455  N5=P1-P(0,1) 

460  FOR  J=1  TO  N2 

481  PRINT  P(0,J)+N5,P(1 ,J)+N5,P<2,JUN5,P(3,J)+N5,P(4,J)*N5 

482  IF  Jsl  THEN  510 

483  REM  Now  we  are  only  looking  at  the  first  and  dominate  harmonic 
500  NEXT  J 

510  PRINT  'Wish  to  store  data  on  DKO : ' ;  \  LINPUT  S$ 

515  IF  S$s ' Y •  THEN  CHAIN  'NSTORE.BAS* 

520  STOP  \  END 


» 


FOUT.BAS 


100  REM  program  FCOMP.BAS 
110  DIM  S%(600),F%(600) 

120  COMMON  M%(3000),A1(4),G(5,7),P{5,7),N,NZ,M1,M2 

130  PRINT  'Wiah  to  compare  the  reaulta  with  original  file’;  \  UNPUT  CS 

135  IF  C$»'N'  THEN  540 

136  DISPLAY  CLEAR 

140  PRINT  'Enter  the  channel  ♦  (or  compariaont';  \  INPUT  M7 
160  FOR  T»M1  TO  M2  \F%(T)»0  \  NEXT  T 
180  FOR  T*M1  TO  M2 
200  FOR  J=1  TO  N2 

220  F%(T)»F%(T)+G(M7,J)»SIN(2*PI*T*J/N+P(M7>J)) 

240  NEXT  J 

260  F%(T)«F%<T)+A1(M7) 

280  S%(T)=M%(T*S+M7) 

300  PRINT  S%(T),F%(T) 

320  NEXT  T 

340  DISPLAY  CLEAR 

360  REGION  OUPPER',1)  \  REGIONCLOWER',2) 

380  GRAPH(,„S%(),2„1)  \  GRAPH!, „F%(),2„2) 

400  LABEL(,'SPIROMETRIC  RESULT  one  breath','.',!) 

420  LABEL(,'FOURIER  ANALYSIS  REPRESENTATION’,'.'^) 

440  PRINT  'The  number  of  harmonica  found  is:',N2 

460  PRINT  PRINT  'CLEAR  THE  SCREEN?';  \  UNPUT  Q* 

480  PRINT  'AGAIN';  \  UNPUT  A* 

500  IF  AJ='Y'  THEN  136 

520  IF  Q$-’Y'  THEN  DISPLAY  CLEAR 

540  CHAIN  'FNORM.BAS* 
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NSTORE.BAS 


100  REM  program  NSTORE.BAS 

120  COMMON  M5l(3000)t  A 1  (4),G(5,7),P(5,7),NtN2,M  1  ,M2 
140  PRINT  'Eater  the  filename  N (tile). DAT;  \  LINPUT  F$ 

160  OPEN  FS&'.DAT*  FOR  OUTPUT  AS  FILE  #5 

180  PRINT  FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA  *•*' 

200  PRINT  #5," 

220  PRINT  #5,'  NORMALIZED  AMPUTUDES' 

240  PRINT  #5," 

260  PRINT  #5, 'SPIROMETER  MAGI  MAG2  MAG3  MAG4' 

270  G1»G(0,1) 

280  FOR  J-l  TO  N2 

300  REM  PRINT  #5,G(0,J)/G(0,1),G(1,J)/G(1,1)>G(2,J)/G(2,1)1G(3,J)/G(3,1),G(4,J)/ 

G(4,l) 

310  PRINT  #5,G(0,J)/G1,G(1,J)/G1,G(2,J)/G1,G(3,J)/G1,G(4,J)/G1 
320  NEXT  J 
340  PRINT  #5," 

360  PRINT  #5/  NORMALIZED  PHASE  ANGLES' 

380  PRINT  #5," 

400  PRINT  #5, 'SPIROMETER  MAGI  MAG2  MAG3  MAG4' 

410  P1*P(0,1) 

415  N5-PI/2-P1 
420  FOR  J*1  TO  N2 

440  PRINT  #5,P(0,J)+N5,P(1  ,J)+N  5,P(2, J)+N  5,P(3,J)+N 5,P(4,  J)+N  5 

450  REM  PRINT  #5,P(0,J|/P1,P(1,J)/P1,P«,.J)/P1,P(3,J)/P1,P(4,J)/P1 

455  IF  J«1  THEN  480 

460  NEXT  J 

480  CLOSE  #5 

500  STOP  \  END 
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Figure  A.l  Schematic  wave  form  illustrating  the  fileplace  -  data 
representation  of  a  typical  forty-five  second  trace. 

By  specifying  various  maxima  and  minima  values,  we  can 
analyze  any  breath(s). 


M  KK*  240  . 
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Appendix  2 

Correlation  coefficients  averaged  over  the  entire  test  duration  for 
each  subject  (four  tables) 


I 


« 


« 


AP  QUIET  BREATHING 


SUBJECT 

S-M] 

S-M9 

S-Mi 

S~Ma 

1 

.9756 

.7029 

.2859 

.7602 

.9604 

.6737 

.4095 

.5166 

2 

.9641 

.0321 

.7187 

.8617 

.9733 

.114 

.8701 

.9320 

3 

.9684 

.5499 

.8772 

.9402 

.978 

.5955 

.9184 

.9685 

4 

.9778 

.6848 

.8313 

.8529 

.9942 

.8666 

.9516 

.8994 

5 

.9287 

.9487 

.7817 

.8928 

.9258 

.8915 

.8529 

.8463 

6 

.9932 

.1206 

.7494 

.7081 

.9923 

-.1543 

.7080 

.7669 

7 

.9575 

-.3899 

.7365 

.7554 

.9795 

-.452 

.7887 

.7360 

8 

.9433 

.4992 

.7954 

.5023 

.9611 

.7998 

.8759 

.6865 

9 

.8081 

.5431 

.9005 

.9243 

.8123 

.8671 

.9632 

.9566 

10 

.9817 

.7019 

.8107 

.9496 

.98 

.6327 

.9501 

.8309 

11 

.9549 

.3118 

.5573 

.8089 

.9000 

.6635 

.491 

.5013 
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AP  FORCED  BREATHING 


SUBJECT 

S-Mt 

S-M? 

S-Mi 

S-M^ 

1 

.9849 

.8845 

.7125 

.6538 

.9825 

.4154 

.5152 

.7858 

2 

.9847 

.3966 

.8795 

.9500 

.9742 

.5829 

.9202 

.8796 

3 

.9562 

.3586 

.9099 

.6506 

.9817 

.0644 

.7709 

.6267 

4 

.9885 

.8884 

.8784 

.6715 

.9873 

.7373 

.8222 

.7403 

5 

.9344 

.9608 

.8360 

.8589 

.9400 

.9646 

.7994 

.8065 

6 

.9924 

.0985 

.6358 

.3156 

.9493 

.1685 

.7157 

.4737 

7 

.9747 

-.1161 

.7514 

.7900 

.9815 

-.4587 

.8371 

.8058 

8 

.9542 

.3257 

.7513 

.5126 

.9672 

.7193 

.9250 

.4977 

9 

.894 

.083 

.8755 

.7737 

.9621 

.0929 

.9284 

.8048 

10 

.9815 

.6555 

.8428 

.6209 

.9914 

.6953 

.8102 

.5193 

11 

.9433 

.4500 

.5396 

.5653 

.9788 

.0986 

.8221 

.8171 

t 
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LATERAL  FORCED  BREATHING 


SUBJECT 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 


S-Ms 

S-M? 

S-Mfl 

.9498 

.8243 

.4687 

.8271 

.9232 

.8563 

-.2496 

.8062 

-.5975 

.5453 

.8145 

.7288 

-.7195 

.6283 

.7250 

.7511 

.7910 

.7788 

.9334 

.7523 

.8450 

.8362 

.9157 

.8708 

.9453 

.8342 

-.7874 

.9471 

.7000 

.8479 

-.6118 

.8523 

.4483 

.3807 

.8277 

.694 

.8778 

.4090 

.3967 

.7428 

.9745 

-.0956 

.813 

.5289 

.9128 

.5022 

.8911 

.5499 

.8807 

.4291 

.9126 

.6256 

.1717 

.8939 

.7823 

.8789 

-.0776 

.8987 

.7360 

.4747 

.0314 

.3034 

.3671 

.6653 

.6533 

-.1091 

.6391 

.7598 

.9323 

-.8566 

.3128 

.7848 

.7259 

-.4663 

.1383 

-.2513 

.5952 

-.6212 

-.7816 

-.9421 

.8375 

-.8987 

-.7784 

-.8823 
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LATERAL  QUIET  BREATHING 


SUBJECT 

S-Ms 

SlMd 

S  m7 

1 

.9260 

.1024 

.0612 

.9182 

.1963 

.2576 

2 

-.7599 

.0007 

.8106 

-.9163 

.1814 

.8882 

3 

.8213 

.4033 

.9653 

.8459 

.3365 

.9423 

4 

.8483 

.8142 

.4001 

.8420 

.7663 

.2381 

5 

.9365 

.3881 

.9498 

.9325 

.8679 

.8903 

6 

.9671 

-.4681 

.8297 

.9563 

-.2238 

.8528 

7 

.7661 

.5974 

.8628 

.6750 

.5631 

.7379 

8 

.1779 

.6302 

.2795 

.514 

.8236 

.3090 

9 

.5011 

-.5108 

.9380 

.659 

-.2317 

.6926 

10 

.9478 

-.7823 

.7529 

.7090 

-.8433 

.6728 

11 

.7596 

.3045 

-.1464 

.3086 

-.8804 

-.4177 

.1375 

.1591 

.5375 

.479 

.9443 

.9572 

.9183 

.9458 

.9254 

.9508 

.015 

-.002 

.3024 

.4607 

.2273 

.0596 

.9391 

.9024 

.8037 

.7586 

-.6463 

-.8461 
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Appendix  3 


The  fourier  analysis  results,  one  typical  breath  for  each  subject. 


NAMY2.DAT 


•••  FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA  ••• 


NORMALIZED  AMPLITUDES 


SPIROMETER 

1 

.0818772 
.0366824 
.0285365 
.01 17908 


MAGI 

1.53173 

.135557 

. 127475 

.0706807 

.0161791 


MAG2  MAG 3 

.277958  1.05395 

.153484  .438685 

.0547998  .190196 

8.14475E-03  .0542658 

7.58024E-03  .0500434 


MAG4 

.243016 

.20207 

.151556 

.116529 

.0901269 


NORMALIZED  PHASE  ANGLES 


SPIROMETER  MAGI  MAG2 

1.5708  1.55904*  1.3449 


MAG3  MAG4 

1.57636  1.87677 


NCGK2.DAT 


•••  FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA  ••• 

NORMALIZED  AMPLITUDES 

MAG4 
.886666 
.139206 

.0537833 
.0309087 
9.70657E-03 

NORMALIZED  PHASE  ANGLES 


SPIROMETER 

1 

.120703 

9.54317E-03 

.012031 

6.71657E-03 


MAGI 

.745678 

.068604 

.0201517 

.0231332 

.0191849 


MAG2 

.0797519 

.0306968 

.0705527 

.0289273 

.023747 


MAG3 

.659867 

.133125 

.  104838 

.0291821 

.0108346 


SPIROMETER  MAGI  MAG2 

1.5708  1.54219  2.4549 


MAG3  MAG4 

1.60014  1.59252 
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NDAT2.DAT 


•••  FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA  ••• 


NORMALIZED  AMPLITUDES 


SPIROMETER 

t 

.  158674 
.0229658 
.0349847 
.0174095 


MAGI 

.676668 

.09221 15 

.0191326 

.0147163 

.0228745 


MAG2 

.202765 

.281772 

.0473572 

.0413403 

.0156769 


MAG3 

.889067 

.276767 

.0825194 

.0667264 

.0310486 


MAG4 

.551774 

.0951901 

6.28071E-03 

.0210197 

.0209508 


NORMALIZED  PHASE  ANGLES 

SPIROMETER  MAGI  MAG2 

1.5708  1.61464  1.68415 


MAG3 

1.5423 


MAG4 

1.43012 


NJAG2.DAT 


FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA 


NORMALIZED  AMPLITUDES 


SPIROMETER  MAGI  MAG2 

1  1.07428  .450423 

.18042  .210806  .210626 

.0598286  .0254997  .0425669 

.0575806  .0417476  .0816732 

.0163215  .023997  .0103082 


MAG3 

1  .02361 

.320047 

.109846 

.135963 

.0301155 


MAG4 

.  190903 

.0730519 

.0476739 

.0318352 

.0224769 


NORMALIZED  PHASE  ANGLES 


SPIROMETER  MAGI 

1.5708  1.53719 


MAG2 

1.820j5 


MAG  3 
1.67781 


MAG  4 
1.63941 
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NJDH2.DAT 


•••  FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA  ••• 

NORMALIZED  AMPLITUDES 

SPIROMETER  MAGI  MAG2  MAG3  MAG4 

1  .308159  .368654  1.28451  1.1958 

.0590098  .027976  .106163  .259992  .129805 

.049087  .0307273  .0438816  .126412  .0659156 

.0124552  .0582475  .0301776  7.07686E-03  .0467637 

.0188293  3 . 32966E-03  3-06696E-03  .0193889  3-98701E-03 

NORMALIZED  PHASE  ANGLES 


SPIROMETER  MAGI  MAG2  MAG 3  MAG4 

1.5708  1.18552  1.60638  1.72337  1.68595 


NKGM2.DAT 


***  FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA  ••• 

NORMALIZED  AMPLITUDES 

SPIROMETER  MAGI  MAG2  MAG 3  MAG4 

1  1.7058  .112815  .657055  .276472 

.176537  .275889  .0844583  .144678  .0588457 

.0380045  .0994377  .0753942  .0952208  .0317776 

.0215485  .0150283  .036331  .0369434  .0122785 

5.48140E-03  .0134287  .035262  .0299376  .0203154 

NORMALIZED  PHASE  ANGLES 

SPIROMETER  MAGI  MAG2  MAG 3  MAG4 

1.5708  1.62188  3.40169  1.67889  1.76207 


NRSH2.DAT 


*»»  FOURIER 

ANALYSIS  OF 

MAGNETOMETER  DATA 

NORMALIZED  AMPLITUDES 

SPIROMETER 

MAGI 

MAG  2 

MAG3 

MAG4 

1 

1.06612 

.908002 

1.90461 

.984701 

.320819 

.310658 

.386747 

.569066 

.205221 

. 109077 

. 125224 

.036669 

.178088 

.0861792 

.Q4934Q7 

.0398869 

.144349 

.092692 

.0625694 

2 .891 36E-03 

.0431166 

.0213328 

.0427264 

.0186506 

normalized  phase 

ANGLES 

spirometer 

MAGI 

MAG2 

MAG3 

MAG4 

1.5708 

1 .63064 

1.2926 

1  .58333 

1 .40485 

NRT02 

.DAT 

»•*  FOURIER 

ANALYSIS  OF  MAGNETOMETER  DATA 

NORMALIZED  AMPLITUDES 

SPIROMETER 

MAGI 

MAG  2 

MAG  3 

MAG4 

1 

.432942 

.537349 

1  .03694 

.444539 

.0970627 

.0260461 

.149473 

.218046 

. 106696 

.038051 

.0315812 

.0828117 

.0965678 

.0354931 

.0218677 

6.65992E-03 

.0321194 

.0221 177 

.0294254 

.0129702 

.0140204 

9.25482E-03 

.0359399 

.042443 

NORMALIZED  PHASE  ANGLES 

SPIROMETER 

MAGI 

MAG  2 

MAG3 

MAG4 

1.5708 

1.68511 

1.8608 

1.65996 

1 .74555 
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NSMV2.DAT 


FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA  ••• 

NORMALIZED  AMPLITUDES 

SPIROMETER  MAGI 

1  .443144 

.111953  .119025 

.0871499  .0470324 

.0586643  .0484116 

•0339576  .0344445 

NORMALIZED  PHASE  ANGLES 

SPIROMETER  MAGI  MAG2  -  MAG  3  MAG4 

1.5708  1.26516  2.04038  1.82385  1.6121 


MAG2 
1.27092 
•  455513 
.26756 
.149006 
.0944689 


MAG  3 
1 .40955 
.30066 
.198715 
.122612 
.0658373 


MAG4 

.835117 

.162786 

. 102874 

.0472155 

.0283885 


NTDK2.DAT 


FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA  ••• 

NORMALIZED  AMPLITUDES 

MAG4 
.450589 
.0604822 
.0511678 
.0475755 
.0418472 

NORMALIZED  PHASE  ANGLES 


SPIROMETER  MAGI  MAG2  MAG3 

1  .850661  .247883  1.22582 

.188684  .213981  .0989264  .353044 

.0458051  .0573182  .0344821  .0998746 

.0660073  .0246433  .0493113  .114086 

9.45356E-03  .0349985  .0326862  .0718003 


SPIROMETER 
1 .5708 


MAGI 

1.85694 


MAG2 
1 .21839 


MAG3 

1.48776 


MAG  4 
.72953 


NWJD2.DAT 


•••  FOURIER  ANALYSIS  OF  MAGNETOMETER  DATA  ••• 
NORMALIZED  AMPLITUDES 


SPIROMETER 

MAGI 

MAG2 

MAG3 

MAG4 

1 

.640421 

2.03128 

1 .66554 

1.08363 

.0666984 

.132997 

.254743 

.279774 

.11212 

.0428129 

.0299331 

.283275 

.0197492 

. 182894 

3.77183E-03 

.0480178 

.17451 

.  176451 

. 149901 

.0205931 

4.52137E-03 

.082091 

.0315118 

.0437652 

NORMALIZED  PHASE  ANGLES 


« 


SPIROMETER 
1 .5708 


MAGI 

1.02658 


MAG2 

2.00377 


MAG3 

1.8433 


MAG  4 

1 .13621 


DATE 

FILMED 


DTIC 


